Non-linear electric transport in graphene: quantum quench dynamics and the 

Schwinger mechanism 



o 
o 

or 

C/3 



Balazs D6ra 1: Q and Roderich Moessner 1 

1 Max- Planck- Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Germany 

(Dated: September 14, 2009) 

We present a unified view of electric transport in undoped clean graphene for finite electric field. 
The weak field results agree with the Kubo approach. For strong electric field, the current increases 
non- linearly with the electric field as E 3 ^ 2 . As the Dirac point is moved around in reciprocal space 
by the field, excited states are generated, in a way analogous to the generation of defects in a 
quench through a quantum critical point. These results are also analyzed in terms of Schwinger's 
pair production and Landau-Zener tunneling. An experiment for cold atoms in optical lattices is 
proposed to test these ideas. 
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The discovery of graphene, a single sheet of carbon 
atoms in a honeycomb lattice (HCL) has triggered intense 
research recently 0, HJ not only because of its potential 
application in future electronic devices, but also because 
of its fundamental physical properties: its quasiparti- 
cles are governed by the two-dimensional Dirac equation, 
and exhibit a variety of compelling (pseudo)relativistic 
phenomena such as the unconventional quantum Hall 
effect 0], a (possibly universal) minimal conductivity at 
vanishing carrier concentration^, Klein tunneling in p-n 
j unctions 0, 0| and ZitterbewegungQ. 

Quantum transport and non-linear responses driven by 
finite external fields represent a genuine non-equilibrium 
phenomenon, giving rise to e.g. dielectric breakdown or 
Bloch oscillations [8(- The quantum aspect of these ef- 
fects is particularly pronounced in reduced dimensions. 
Therefore, two-dimensional Dirac electrons in finite elec- 
tric fields, the subject of this work, provide a fascinating 
setting for studying these issues. 

A simple picture of electronic transport in a finite elec- 
tric field is drift transport as considered by Drude: car- 
riers move ballistically (p = eEt) until they change their 
momentum by a scattering process, replacing the time 
t by the appropriate scattering time. The special fea- 
tures of Dirac electrons relevant for transport in finite 
field include: (i) their velocity is pinned to the "light 
cone" Fermi velocity, vp, (ii) relativistic particles un- 
dergo pair production in strong electric fields, as pre- 
dicted by Schwinger and (iii) a uniform electric field 
modifies locally the geometry of the Fermi surface by 
moving the Dirac point around in momentum space (Eq. 

Since massless Dirac electrons can be thought of as 
being critical, this can lead to the production of excited 
states, and should leave its fingerprints on transport in 
finite electric fields. 

The Landau-Zener (LZ) dynamics, describing the 
(avoided) level crossing in a two level system [loj. repre- 
sents the natural language to discuss Klein tunneling 6, 
HH in graphene, and is ultimately connected to defect for- 
mation and quench dynamics through quantum critical 



points [13], described by the Kibble-Zurek mechanism [13, 
1141 ] of non-equilibrium phase transitions. Applying these 
ideas to graphene allows us to analyze the real time dy- 
namics of the current, after switching on the electric field, 
and to identify the various crossovers (summarized in Ta- 
bic UJ). Electric transport depends sensitively on the fre- 
quency, temperature, electric field and scattering rate (ui, 
T, E, r), and the obtained current depends strongly on 
how the (lo,T,E,T) -> limit is takenjj [lj|. Our re- 
sults follow from taking the to = T = limits in a finite 
electric field. 



classical 
t < h/W 


Kubo 


Schwinger/Kibble-Zurek 


h/W < t < y/h/vpeE 


s/h/vpeE Ci< t B ioch 


jm ~ Et 


3*~E 


j x ~ tE 3 ' 2 



TABLE I: Temporal evolution of the non-equilibrium cur- 
rent for clean graphene. Bloch oscillations show up for 
t > tBioch ~ h/eaE\e§ with a the HCL constant. 



We focus on the 2+1 dimensional Dirac equation in 
a uniform, constant electric field in the x direction, 
switched on at t — 0, through a time dependent vector 
potential as A(t) = (A(t), 0,0) with A(t) = EtQ(t). The 
resulting time dependent Dirac equation, describing low 
energy excitations around the K point in the Brillouin 
zone for clean graphene, is written as 



H = v F [a x (p x - eA(t)) + a y p v ], 

ihd t %(t) = H%(t), 



(1) 



where vp ~ 10 6 m/s is the Fermi velocity of grap- 
hene, and the Pauli matrices (a) arise from the two 
sublatticesQ of the HCL. Due to this (pseudo)spin struc- 
ture, Eq. {TJ represents a natural platform to study LZ 
dynamics as well. It is convenient to perform a time 
dependent unitary transformation first [lj|, which diago- 
nalizes H , and brings us to the adiabatic basis in the LZ 
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language [lfj as 



U 



exp(—ip/2) 
exp(itp/2) 



exp(— itp/2) 
— exp(itp/2) 



U+HU = a z e p (t), where 



e p (t) = v F J(p x - eA(t)) 2 + p 



(2) 



(3) 



and tan tp — p y j (p x — eA(t)). With this, the transformed 
time dependent Dirac equation is given by 



iWt%{t) = 



hvpp y eE 



Hi*) 



^(t = 0) = (0,l), 



(4) 
(5) 



where &p(t) — U$> p (t), and the off diagonal terms in the 
Hamiltonian arise due to the explicit time dependence of 
the unitary transformation (—iU + dtU), and the initial 
condition corresponds to zero temperature and half fill- 
ing. The main advantage of the unitary transformation is 
that the resulting equation clearly distinguishes between 
positive and negative energy states. 

The current operator in the original basis is ob- 
tained through the equation of motion as j x = —ev F a x . 
After the unitary transformation, it reads as j x = 
— ev F (a z cos tp + a y sirup). In the presence of the elec- 
tric field, the expectation value of the current is finite. 
By denoting = (a p (t), f] p (t)), 



(jx) P (t) = -ev F [cosp(\a p (t)\ 2 - \(3 p (t)\ 2 )- 
+2smtpRe(ia p (t)p;(t))] . 



(6) 



The first term is the current from particles residing on 
the upper or lower Dirac cone, while the second one de- 
scribes interference between them, and is responsible for 
Zitterbewegung. Using QED terminology, the first and 
second term is referred to as conduction and polariza- 
tion current, respectively (l7|. In condensed matter, these 
are called intraband and interband contributions, respec- 
tively. Due to charge conservation, |a p (t)| 2 — \/3 p {t)\ 2 = 
2\a p (t)\ 2 — 1. The interference correction also simplifies 
since 



d t \a p (t)\ 2 =2Re(a p {t)d t a* p (t)). 
By using the transformed Hamiltonian, Eq. Q 

. hvpp y eE 



(7) 



hd t a*(t) = ie p (t)a*(t) - i- 



/$(*)> (8) 



consequently 



Vpp y eE 



ft|ap(i)r = - Re(ia p (t)p;(t)), 
e p\ l ) 



(9) 



since Re(i|a p (i)| 2 ) = 0. As a result, the expecta- 
tion value of the current only requires the knowledge of 
n p (t) = \a p (t)\ 2 as 



(jx) P (t) = ~ev F 



vf{Px - eEt) 



(2n p (t) - 1)- 



-2^-d t n p (t) 
VFeE 



(10) 



The term independent of n p (t), namely evp(p x — 
eEt) I e p (t) vanishes at half filling after momentum inte- 
gration. In QED, this originates from charge conjugation 
symmetry [17} . while in graphene, it is obtained by taking 
the full honeycomb lattice into account as in Ref. Hit 

For t < 0, the upper/lower Dirac cone is empty/fully 
occupied. The quantity n p (t) measures the number of 
particles created by the electric field in the upper cone 
through Schwinger's pair production In graphene, 
instead of particle-antiparticle pairs, electron-hole pairs 
are created. Therefore, the basic quantity to determine 
transport through graphene is n p (t). We start by ana- 
lyzing its behaviour at weak electric fields perturbatively. 
In this case, we can set E = in Eq. ((4|) except in the 
numerator of the off-diagonal terms, and obtain 



,(*) = 



(eEhpy)'- 
4«||p| 6 



■ sin 



/ v F \p\t 

v ft 



(11) 



which is valid except in the close vicinity of the Dirac 
point (i.e. \p\ 3> eEt), and \p\ 
this into Eq. dTT 



p\. Plugging 

the first term is already second order 
in the electric field, and does not contribute to linear 
response. The second (polarization) term gives, taking 
valley and spin degeneracies into account 



e 2 E 

{jx) = 2^n 



dp 



sin(2v F pt/H) 



e 2 

Tn E > 



(12) 



and the dc conductivity is a = j/E = e 2 it/2h, in 
accordance with Ref. [l8(. This is the value of the 
ac conductivity at finite frequencies obtained from the 
Kubo formula]!, 0| and measured also [ijj], and since 
the model does not contain any additional energy scale, 
which would change the value of the ac response down 
to u — > 0, the same value for the dc conductivity sounds 
plausible. Within our approach, the small field response 
is dominated by Zitterbewegung corrections. The ultra- 
short time transient response (tW <§; h with W the band- 
width) is fully classical. Expanding Eq. (fTTj) . we obtain 



p 2 

(j x ) P (t)=e 2 v F j^Et, 



(13) 



independent of H. The current rises linearly with time 
after the switch on as (j x )(t) = 4e 2 EWt/h 2 . The very 
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same result follows from a classical Hamiltonian, H c i 
■p 2 . The Hamilton equation is 



VfJ (Px - eEt) 2 



d t x 



dH c , 
dp x 



,(p x - eEt) 

Hrl 



(14) 



which gives for the classical current, j c i(p, t) — —edtX, at 
short times as in Eq. (fT5|) . Dirac particles can therefore 
be accelerated as dfx = eE/m xx at short times, in accord 
with Newton's equation after defining their effective mass 
as l/m xx = d 2 H/dp 2 x = v F p\l\p z \- 

For the general time and electric field dependence, Eq. 
(HJ) can be solved analytically [l7l l2(ij using the parabolic 
cylinder functions, which do not immediately yield a 
transparent analytical expression for the electric current 
for arbitrary electric field and time. To investigate the 
strong field, long time (specified in Eq. fl8|)) response 
of Dirac electrons, we use the asymptotic expansion of 
these eigenfunctions[l6|, 03, or equivalently we can rely 
on the WKB approach [2 lj to determine n p (t) through 
the barrier penetration factor, similarly to narrow gap 
semiconductors [2^ . As a result, we get 



n p (t) = Q(p x )Q(eEt - p x ) exp 



heE 



(15) 



which is the celebrated pair production rate by 
Schwingerj^, [13], a manifestation of Klein tunneling^, 
and also the LZ transition probability [lOj between the 
initial and final levels. More precisely, the conditions for 
applicability are (p Xl eEt — p x ) ^> \p y \. This expression 
can be transparently understood invoking LZ physics. 
Two levels at ±p x , weakly coupled by p y level cross with 
time, ending up at ±(p x — eEt). The transition is com- 
pleted when both the initial and final levels are well sep- 
arated, in which case the mixing between them is given 
by Eq. ([IB}, as plotted in Fig. □ 

Putting Eq. (fTS"|) into Eq. fllOjl . the current is domi- 
nated by the conduction (intraband part) as 



<?«>(*) 



2e 2 E v F eEt 2 



(16) 



which increases linearly with time, similarly to normal 
electrons in a parabolic band. However, the origin of the 
time dependence is completely different: it stems from 
the increasing number of pairs due to pair production a la 
Schwinger, each contributing with the same velocity vf, 
as opposed to the continuously accelerated fixed number 
of normal electrons in strong fields. 

The total number of particles and holes created is 



N(t) 



j^Jdpn p (t) = - 



2eE 



Vp 



eEt 2 



(17) 



which leads to Eq. (flUj) via (j x )(t) — ev F N(t), and is re- 
lated to the quench dynamics through a quantum critical 



point (QCP)[23[ as follows: Eq. JTJ) can be diagonalized 
at every instant with eigenenergies in Eq. ([3]): the Dirac 
point moves continuously in momentum space with loca- 
tion p = (eEt, 0), which results in defect (excitation) pro- 
duction. The spectra from Eq. ([3]) can be considered as 
an ensemble of 1+1 dimensional initially gapped systems 
(labeled by p x ) driven through a QCP. The initial energy 
gap is given by wfIpxI, the one dimensional momentum is 
p y , and the quench is applied as vp(p x ~eEt). For a given 
p x , during the temporal evolution, the gap vanishes at 
the instant t — p x /eE, which defines the QCP, and reap- 
pears with increasing time. The dynamics close to the 
QCP (characterized by d = z = v = l[23j]) is necessar- 
ily non-adiabatic (impulse) due to the divergence of the 
relaxation time and the finite quench time ~ l/ej^fl^]. 
The Kibble-Zurek mechanism 131 1 141) predicts a scaling 
form for the defect formation asQ J5<W(*"+l) = 
for a given 1+1 dimensional system. 

However, defect production occurs only upon complete 
non-adiabatic passage through the QCP. At a given time 
t, this holds for <C p x <C eEt, so the number of 
quenched systems scales ~ tE. Combining these, the 
Kibble-Zurek mechanism thus also predicts the tE 3 / 2 
scaling of the total defect density for Eq. ([1]), similarly 
to Eq. (|17| , linking the non-linear transport in graphene 
to critical phenomena. (Quantum critical transport from 
a different perspective was already studied in Ref. I25I ) 



+ eEt. 
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FIG. 1: (Color online) Left panel: visualization of the tempo- 
ral evolution of the LZ dynamics. Right panel: schematic pic- 
ture of the current-electric field characteristics for graphene. 
Interband transitions are overwhelmed by intraband ones with 
increasing electric field, and the character of the measured 
current changes from polarization to conduction contribution. 



Therefore, the low field, perturbative response is dom- 
inated by interband contributions, and can be regarded 
as a manifestation of Zitterbewegung. With increasing 
field, a large number of electron-hole pairs are created, 
and intraband processes take over, producing non-linear 
transport. The crossover is determined by the dimcn- 
sionless time-scale, after comparing our system to the 
LZ model asllOl 



v F eEt 2 



(18) 
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For T cross < 1, no level crossing occurs, and we can use 
perturbation theory to estimate the current, therefore we 
are in the Kubo regime. The Kibble-Zurek mechanism 



defines the freeze-out time|12l. 1 141 ] by the instant t when 
the system leaves the adiabatic regime and enters into 
the impulse one, namely i — H/vFeEt, and the Kibble- 
Zurek form of the defect density requires complete transit 
through the QCP, T cross > 1). In the LZ 

language, level crossing is completed for T cross ^> 1, the 
number of pairs created is non-perturbative in the electric 
field, and we can use the probability of LZ tunneling for 
the current. 

So far we have discussed the real time evolution of 
the current after the switch-on of the electric field, sum- 
marized in Table |U In ideal clean graphene, for long 
enough times, Bloch oscillation would set in due to the 
underlying HCL structure. In reality, the time t must 
be replaced, in the spirit of the Drude theory, by an ap- 
propriate scattering time [HI (due to phonons or impu- 
rities), t sc , or in ballistic samples, by the ballistic flight 



time from the finite flake size, Tf, = -^/^fIIII- The ob- 
servation of non-linear electric transport requires, from 
Eq. (p~8|) . an electric field as 



E> E c = H/v F eT 2 



(19) 



where r = min(Ts C , r^, ta) is the shortest of the addi- 
tional restricting time scales (with ta defined below). 
Ballistic transport on the (sub)/im scale implies r ~ 
0.1 - 1 ps, giving E c ~ 10 3 - 10 5 V/m0. The mea- 
sured current is expected to show a change of slope as a 
function of the electric field in the crossover region, and 
an extended electric field window would be required to 
reveal the non-integer exponent, as shown in Fig. [TJ It is 
important to emphasize that in both regions, the current 
is related to n p (t), thus even the linear response regime 
witnesses pair production. 

In the presence of a small mass gap, the above re- 
sults need to be modified. The perturbative regime is 
characterized by exponentially activated behaviour due 
to the gap, and the current is exponentially suppressed 
at low temperatures (T <C A) as j ~ Eexp(— A/T), as 
in normal semiconductors. On the other hand, for strong 
electric field, we can still use the analogy to LZ tunneling 
as 



<J*>(<) 



2e 2 E v F eEt 2 



■ exp 



irA 2 



hvFeE ) 



(20) 



Non-linear transport sets in for E > it A 2 /Hvfg, which 
defines a new timescale for E c as ta = h/ Ay/n. 

In general, the non-linear current for d + 1 dimen- 
sional (d — 1, 2, 3) Dirac electrons^ is (j x )(t) ~ 
tE (d+i)/2 ex p(— ttA 2 /%v F eE). For d = 1, a good real- 
ization would be carbon nanotubes (rolled up graphene 
sheet), whose "non-linear" response is still linear (j ~ E), 
only the non-trivial exponential factor with a possi- 
ble gap reports about non-perturbative effects [27]. The 



d — 3 case could be realized among the bulk electrons of 
Bi, possessing a band-gap ~ 0.015 eV. 

These results are also relevant for other systems with 
possible Dirac fermions such as the organic conductor [28j| 
a-(BEDT-TTF) 2 I 3 with a tilted Dirac cone. Dirac 
fermions can be realized in cold atoms in an appropriate 
optical lattice (half filled HCL, Kagome and triangular 
lattices), without any source of dissipation or scattering. 
The momentum distribution, Eq. (|15|) reveals the effect 
of the driving electric field before Bloch oscillations set 
in0]. The pairs created increase the energy of the system 
as ~ t 2 E 5 / 2 , which, together with the momentum distri- 
bution of Bloch states, can be measured after releasing 
the trap. This could be a first direct experimental obser- 
vation of the Schwinger mechanism as well. 

Note added Recently we became aware of a related 
work [29j|. Overlapping results are in agreement. 
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